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Abstract 

Numerical simulations of the evolution of a global topological defect field have two char- 
acteristic length scales — one macrophysical, of order the field correlation length, and the 
other microphysical, of order the field width. The situation currently of most interest to 
particle cosmologists involves the behaviour of a GUT-scale defect field at the epoch of 
decoupling, where the ratio of these scales is typically of order 10 50 . Such a ratio is unre- 
alisable in numerical work, and we consider the approximations which may be employed 
to deal with this. Focusing on the case of global texture we outline the implementation 
of the associated algorithms, and in particular note the subtleties involved in handling 
texture unwinding events. Comparing the results in each approach then establishes that, 
subject to certain constraints on the minimum grid resolution, the methods described are 
both robust and consistent with one another. 



PACS numbers: 98.80.Cq, 98.70.Vc 



1 Introduction 



Topological defects as a class of phenomena naturally emerge in the very early universe from 
the symmetry breaking phase transitions associated with Grand Unifying extensions to the 
standard SU(3)c <8> SU(2)l®U(1)y model of particle physics jjj. A minimal requirement for 
any theory incorporating such defects is that the energy density in the defect field does not 
dominate the universe today. In certain cases the defect field is believed to evolve to a scaling 
regime, wherein its energy density remains a fixed small fraction of the total. Such defects can 
then be invoked as a possible source of the primordial density perturbations believed to have 
seeded cosmological structure formation and observed as fluctuations in the cosmic microwave 
background (CMB) by the COBE satellite 0. The essential observational constraints on any 
model of structure formation based on topological defects — most notably invoking either 
gauged cosmic string || or global texture j| — are therefore of the large scale structure of 
the local universe essentially as it is now, at a time t ~ 10 18 s after the big bang, and of the 
CMB radiation from the epoch of matter /radiation decoupling at t ~ 10 13 s. By contrast 
the assumed defect-forming Grand Unified Theory symmetry breaking typically occurs at 
t ~ 10 -35 s. Numerical simulations of the evolution and effects of the defect field have 
therefore come to play an essential role in bridging this gap of some 50 orders of magnitude 
in time. 

A crucial factor determining the type of numerical simulation used is whether the sym- 
metry broken at the phase transition is gauged or global. For a broken gauged symmetry 
any vacuum field configuration which is not topologically stablised can simply be trivialised 
by a suitable gauge transformation, and all that need be considered is the behaviour of the 
defects themselves which can be treated as discrete objects. For a broken global symmetry 
however topologically unstable vacuum field configurations are only trivialised by correlation 
as the field evolves. Since the correlation length of a global defect field is typically taken to 
be of order the horizon scale its evolution must be followed in detail up to this scale. Cos- 
mological global defect field simulations thus contain two fundamental length scales — one 
macrophysical, of order the horizon size, and the other microphysical, of order the inverse 
mass of the field — which we can quantify via their ratio as 

t n~ l 
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For a GUT-scale defect at the epoch of decoupling 1Z ~ 10 50 and our characteristic length 
scales again differ by 50 orders of magnitude. Since it is not yet possible to simulate such 
a range of scales simultaneously we must adopt one of two possible approximation schemes. 
One approach is to reduce the field's mass scale, and hence increase its width, to a level which 
is numerically resolvable || ^, [7|, g|. Since the field's width characterises the size of the defect 
core this may be dubbed the expanded core (XCORE) approximation. The alternative is to 
take the field's mass scale to infinity and hence its width to zero, removing the microphysical 
scale from the problem H, |l^, 11, 12]. Provided this is done self-consistently the field's 



equations of motion then reduce to those of the non-linear sigma model (NLSM). 

In this paper we consider the implementation of numerical simulations of the evolution 
of global texture in each of these approaches. In each case the approximated field mass is 
very different from the physical mass, and it is essential to determine the conditions under 
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which the resulting simulated evolution still reflects the true behaviour of the field. In the 
XCORE case some progress has already been made [15] by finding the minimal grid resolution 
necessary for some of the properties characterising the field's evolution to asymptote. Here, 
however, we are able to go further; since the two approaches involve making diametrically 
opposed approximations consistency between them would be strong evidence that each is 
close to the true field evolution. Our primary focus therefore is on comparing results from 
each simulation approach and testing their mutual consistency. 

As well as being one of the more promising candidate defect models, global texture also 
displays some unique features in its evolution. In particular although the texture field is glob- 
ally topologically trivial, being everywhere on the vacuum manifold with no regions trapped 
in the false-vacuum phase, it nonetheless admits configurations which are in some sense lo- 
cally topologically non-trivial. In more than two spatial dimensions such configurations are 
unstable, with their field gradient energy concentrating (by analogy with Derrick's theorem 
|ji~3|| ) until the local energy density is sufficient to pull the field off the manifold and through 
the false-vacuum, 'unwinding' this locally non-trivial topology. Since such unwinding events, 
where the field leaves the vacuum manifold, are where the different approximation schemes 
are most likely to diverge, they will provide the most stringent test of consistency. 

Having established our notation we commence by detailing the continuum form and discre- 
tised implementation of the equations of motion in each approximation scheme. We include in 
this some discussion of the identification of unwinding events and an analysis of the minimal 
grid resolution necessary for the approximations to model such events accurately. Working 
within these constraints we are then able, at least in flat space, to check the simulations 
against one another. For a number of initial field configurations, both with and without 
unwinding events, we consider the field's point-by-point evolution, its gross features, and the 
CMB anisotropies it generates in each approximation. 



2 Notation And Conventions 

Except where specified otherwise we consider the evolution of a 4-component real scalar global 
field in 3 + 1 dimensional flat space with metric signature (—,+,+,+). This field evolves 
in a potential of the form 

y($) = K(|$[ 2 -$o) 2 (2) 

which is defined such that above some critical temperature < and the potential has an 
0(4) symmetry about its minimum at $ = 0, whilst below it ^ > and the symmetry is 
broken to 0(3) about the new minimum which lies somewhere on the entire S3 = 0(4)/ 0(3) 
manifold of degenerate vacuum states defined by |$| = <3? - This is the simplest symmetry 
breaking to generate global texture. 

In discussing texture configurations it is often useful to parameterise the field at some 
point in space given by the spherical polar co-ordinates (r,8,ijj) at some time t as 

^>(r,9,ip,t) = |3>|(cos X, sin X cos 6, sin X sin 6 cos \I>, sin X sin 6 sin*) (3) 

where the group space angles < ^ < 7r, < < 7r and < ^ < 2n locate the projection 
of on the vacuum manifold. In general X, O and \& are functions of position and time, 
although spherically symmetric field configurations can be realised in a particularly simple 



2 



way by identifying two of the group space angles with the two physical space angles, = 6 
and \P = if). 

For convenience we further normalise the field, <3? — > 3>/3> - All simulations are now 
performed on lattices with grid spacing 5x and using time steps 5t and with the speed of 
light set to c = 1 so that grid distances and times are numerically interchangeable. The 
discretised field is labelled f»j )n with i-indices spatial and n-indices temporal. 

3 Implementation 

3.1 Equations Of Motion 

The full texture Lagrangian 

c = - l -d^d^-v Q m 2 -if (4) 

gives the equations of motion of the field as 

dV 

□* = -- (5) 

In flat space these reduce to 

<t> - V 2 $ = -4U (|$| 2 - 1)$ (6) 



Using a standard staggered leapfrog method [14] these can be discretised to second order to 
yield 

*i,n+l/2 = *i,n-X/2 + (V 2 $i,„ - 4V (*? n - l)* i>n ) 5t 

$i,n+l = *i,n + *i,n+l/2<ft ( 7 ) 

These are the equations used in the XCORE approximation. 

Alternatively, imposing the NLSM approximation and restricting the field to remain on 
the vacuum manifold at all times, the constrained Lagrangian can be written 

c=-\d^d^-^m 2 -i) (8) 

where a is a Lagrange multiplier. The equations of motion now become 

□ $ = cr$ (9) 

giving 

o- =(□$).$ (10) 
Using the constraint to solve this equation for a the flat space equations become 

«£ - v 2 $ = -(6 2 - v$ 2 )$ (11) 

Again discretising to second order these give 

$i >n+ i = 2$ i)H - $i, n _i + V 2 $,, n ^ 2 - (6? - V$ 2 n )5t 2 $i, n (12) 
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Following Pen et al (T^] this can be cast in the form 

*i,n+l = X $i,n + S$i,n (13) 

where 

A = 1 - ($?„ - V$? n )<5t 2 (14) 
Using the constraint |$j iW+ i| = 1 to solve for A now gives 

A = -$*»,„.*<,„ ± y/l - + ($i, n .5$i, n ) 2 (15) 

and we are able to circumvent the numerically troublesome squared derivative terms. Note 
that to first order in 5t equation (|l3|) is 

= ±* iin + O(tffc) (16) 

with the choice of sign coming from the choice of root in equation (|i~5|) ; taking the negative 
root, so that <J?i jn +i ~ — &i.m therefore explicitly introduces an unwinding event at the given 
gridpoint. 

3.2 Unwinding Events 

Evolving the texture field on a discrete lattice there are two types of unwinding event. If 
the initial field conditions are sufficiently symmetric about some gridpoint that the group 
space angles there are time-independent then an unwinding can occur at that gridpoint. 
However for more random initial field conditions the unwinding site will not lie exactly at a 
gridpoint. To investigate the properties of individual texture unwindings, as well as to test 
the accuracy of our simulations generally, we would like to be able to identify both types of 
event unambiguously. 

The full texture field equations of motion include unwinding events dynamically. Reducing 
the mass scale of the field allows it to be pulled off the vacuum manifold more easily, so that 
making the XCORE approximation increases the volume of space leaving the vacuum and 
forming the defect core. The constraint on the field mass is then that it should be sufficiently 
small that any unwinding event will pull the field below some threshold value — typically 
taken as |$i,n| < 0.5 — throughout a volume containing at least one gridpoint. Conversely the 
mass must also be large enough that the field is only pulled below this threshold at unwinding 
events. However, since we are trying to maximise the range of length scales, and hence require 
as large a mass as possible, only the upper limit serves as a serious restriction. A further 
constraint on the field's mass comes from the upper limit on field gradients imposed by the 
finite grid spacing; if the width is too small then even unwindings occuring at gridpoints can 
be inhibited by the resulting upper limit on the field's gradient energy density. A detailed 
investigation of these constraints |l5| leads to a practical lower limit on the field width 

W>Q.2h5x (17) 
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Subject to this, and the grid resolution conditions discussed below, the threshold prescription 
serves to identify unwindings in the XCORE approximation. 

Since the NLSM approximation constrains the field to remain on the vacuum manifold at 
all times it is unable to admit unwinding events in the continuum. However, as demonstrated 
above, the discretised equations of motion include a means by which an unwinding can be 
explicitly re-introduced at a given gridpoint. All that remains is to determine the appropriate 
criteria under which such an unwinding should be incorporated in the field's evolution. For 
diagramatic simplicity in the remainder of this section we consider the case of a 3-component 
NLSM texture field in 2 + 1 dimensions where the vacuum manifold is S^. 

The most obvious approach is to introduce an unwinding at any gridpoint in the immediate 
vicinity of which the field covers more than half of the vacuum manifold — ie. if we consider 
the bisection of the manifold by the closed boundary on it given by joining the nearest 
neighbours to any gridpoint along geodesies then we introduce an unwinding whenever the 
field at the gridpoint lies in the larger of the two sub-sections. To date the only attempt to 
implement this criterion (1(1, 11] involves calculating the mean value of the (here 4) nearest 



neighbours to a gridpoint, $i >n , and introducing an unwinding if and only if <£j jn • < 0. 
This is equivalent to replacing the actual boundary on the manifold with an 'averaged' one 
defined by the intersection of the manifold with the plane perpendicular to <£j jn (figure 1). 
Note that without loss of generality we can always take §i^ n to point 'south', whereupon 
®i,n ' &i,n < if and only if <£j jn and $i >n lie in opposite hemispheres. 
In the symmetric case it is straightforward to show that 

$i,n OC $i, n (18) 

and hence that <J>j )n and &i. n lie in opposite hemispheres if and only if more than half of the 
vacuum manifold is covered locally (figure 2). 

For non-symmetric configurations this approach breaks down in a number of ways: 

• since in general <&i tH gt ^i >n it is simple to construct configurations where <3?j )n and $i >n 
lie in the same hemisphere but which should nonetheless be unwound (figure 3). As is 
clear from this example the correct condition for identifying an unwinding is instead 
($ i>n - § i>n ) • §i. n < 0. 

• if the field at the gridpoint in question lies between the true boundary and the 'average' 
boundary then even the correct condition either mistakenly identifies or fails to identify 
an unwinding configuration (figure 4). 

• since in general the unwinding site is not at a gridpoint even if the algorithm correctly 
identifies a configuration that should unwind what we are doing in effect is unwinding 
the gridpoint closest to the actual unwinding site. However in the NLSM approxima- 
tion only at the unwinding site itself should the field ever leave the vacuum manifold. 
Moreover, even if the NLSM condition is relaxed and a small region of space close to the 
unwinding site (and including the gridpoint) is allowed to leave the manifold it is only 
at the unwinding site itself that the field is approximately inverted, <&i, n +i ~ — *&i,n- 
This approach therefore incorrectly evolves the field at the gridpoint, and can thereby 
generate a repeated unwinding (figure 5). 
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In practice unless the unwinding site is fixed at a gridpoint (when the required symmetry 
avoids the above difficulties) it is unnecessary to explicitly unwind the field at all, since the 
geodesic assumption (that between neighbouring points the field follows a geodesic on the 
manifold) implicitly unwinds any configuration locally covering more than half the mani- 
fold anyway. However the problem of unambiguously identifying unwindings in the NLSM 
approximation in order to study their properties remains unsolved. 



3.3 Grid Resolution 

Given the above microphysical length scales, the resolution afforded by a particular simu- 
lation lattice can be quantified by the number of gridpoints TV; corresponding to the initial 
macrophysical length scale. For comparison between different approaches we need a length 
scale which is well-defined and which can be meaningfully implemented in all simulations 
- for example in flat space the horizon is not well-defined, whilst symmetric initial field 
configurations have an infinite correlation length. 

A necessary condition for any field configuration to unwind that to do so must reduce the 
gradient energy in the configuration, and hence that it should cover more than half of the 
vacuum manifold before unwinding, and by corollary less than half afterwards. Therefore, 
as in previous work 0, §, we adopt as our macrophysical length scale the radius of a sphere 
within which no more than half of the manifold can possibly be covered initially. For singly- 
wound spherically symmetric configurations |?J the manifold covering about the origin is 
given (in the ansatz of equation (|3|)) by 

Q(r) = i (x (r) - ~ sin 2X (r)) (19) 

and our length scale is the radius at which Q = 0.5 initially. For interpolated random 
configurations ||] it is half the distance between gridpoints at which the initial field values 
are randomly assigned. 

Figure 6 illustrates the variation in the topological charge causally associated with the 
unwinding of a spherically symmetric configuration in flat space Q(t uw — t Q ) (where t uw and 
t Q are the unwinding and initial times respectively), with the grid resolution Ni in each ap- 
proximation. To get within 10% of the asymptotic result in both simulations simultaneously 
a lower limit on the number of gridpoints per half-manifold radius in flat space is Ni > 5, 
and to within 5%, Ni > 10. However it should be stressed that these results only apply 
to spherically symmetric configurations in flat space. Experiments using the XCORE ap- 
proximation in expanding backgrounds and with random configurations demonstrate that 



in such circumstances the minimum resolution is increased [15], and wherever possible the 



simulations presented here satisfy the tighter constraint from that work, iVj > 16. 



4 Comparison 

The evolution of an initial field configuration under each of the approximations can be com- 
pared in a number of ways; here we consider three levels of consistency. Firstly we can 
calculate point by point statistics measuring the correlation of the field over the lattice at 
every time step. Alternatively we can consider the consistency of any gross features of the 
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simulations, such as the times and locations of any identifiable unwinding events. Finally 
we can compare observable consequences of the field evolutions, such as their induced CMB 
anisotropies. Clearly these levels of consitency are hierarchical; if the simulations are con- 
sistent point by point then their gross features and physical consequences will be identical. 
Similarly even if their point by point evolution varies, provided the gross features are the 
same then we might expect the induced CMB anisotropies to be so too. 

4.1 Pointwise Statistics 

Given that the vacuum manifold is the 3-sphere |$| = 1 a simple measure of the correlation 
of the field at any point and time is given (in an obvious notation) by 

= $* n • $£ n (20) 

However we know at the outset that differences will emerge as the simulations progress simply 
because the field is able to leave the vacuum manifold in the XCORE but not in the NLSM 
approximation. As a further comparison therefore we can factor out this effect by comparing 
the normalised fields, projecting the XCORE field back onto the vacuum manifold everywhere 

4n = J?, • *?„ (21) 

I i,n I 

In this way we can determine the relative importance in any absence of correlation of differ- 
ences in the field positions' magnitudes and of differences in their orientations. 

Symmetric configurations are artificially self-correlated on all scales at all times, possibly 
enhancing their cross-correlation too, so we focus on the evolution of random configurations 
alone. Using a catalogue of random initial configurations admitting unwinding events || we 
take a sample of 22 configurations, equally divided between ones with and ones without an 
unwinding. These are then simultaneously evolved on a pair of 48 3 lattices for 24 grid time 
units. Running the simulations simultaneously restricts us to half the usual grid size (and 
hence resolution) here, so the results should be taken as upper limits which would be reduced 
were it possible to work at full resolution. We calculate as our correlation statistic the mean 
values of each measure over the simulation lattice at each timestep 

e» = ^E&,» e. = ^E4n (22) 

i=l i=l 

Exact correlation gives a statistic of 1, anti-correlation -1, and no correlation 0. In all 
simulations and at every timestep we find £ n = £ n = 1.0. 

To test of the degree of this consistency we also note the number and distribution of 
gridpoints at which either of the correlation measures falls below 0.9. In simulations of 
both unwinding and non-unwinding configurations gridpoints whose fields are less than 90% 
correlated are found localised in time and space about events during which the field is pulled 
significantly off the manifold in the XCORE simulations. This is true both for the known 
unwinding events, and for the non- unwinding events (at which insufficient gradient energy 
collapses to unwind the field) which may be identified post hoc by the characteristic CMB 
anisotropy pattern they induce Jig ]. However, these imperfectly correlated regions are very 
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small; even the largest region involves less than 0.1%, and a typical region less than 0.01%, 
of the simulation gridpoints. 

Figure 7 shows the variation in the number of gridpoints whose fields are less than 90% 
correlated under each measure as the simulation progresses for a typical configuration, in 
this case including an unwinding at grid time 13.9. We see that for the full field the number 
of imperfectly correlated points peaks just before and after the unwinding event, but that 
only a small fraction of this discrepancy is due to differences in orientation, and that these 
differences peak at the unwinding itself. The sequence of post-unwinding peaks in the full 
field correlation measure is to be expected given the suppression of the damping of radial 
oscillations by the reduction in the field's mass in the XCORE approximation. Note that, 
except for the first, these peaks involve a single point; explicit examination of the simulation 
confirms that this is indeed the site of the unwinding event. 



4.2 Gross Field Features 

As shown above the only even marginally uncorrelated sections of the simulations are asso- 
ciated with unwinding or non-unwinding events. We now consider the extent to which the 
gross features of these events, such as their time and place, reflect this and differ between 
the simulations. The cross-correlation results already limit the possible differences; in each 
case the event must occur within the imperfectly correlated region of simulation spacetime, 
so their locations cannot differ by more than the size of this region. 

As a more precise test we can attempt to identify and thereby locate each simulations' 
unwinding and non-unwinding events explicitly. The unwinding events are identified in the 
XCORE simulation as occuring at the gridpoint at which the magnitude of the field falls 
below some threshold, |<3?| < 0.5, and at the time half way between its first falling below the 
threshold and subsequently exceeding it again. In the NLSM simulations we use the criterion 
discussed above, ($i, n — $»,n) • &i,n < 0, whilst acknowledging its weaknesses. Note that the 
criterion is used simply to identify candidate unwinding events, and the equations of motion 
are unchanged, with the positive root being taken in equation di~5|). Where unwindings are 
successfully identified (in all but one case) they are found to occur at the same gridpoint and 
within 0.5 grid time units of their XCORE counterparts. 

Identification of non-unwinding events is less precise; in both XCORE and NLSM simula- 
tions we look for the characteristic CMB signature — a sharply peaked cold spot followed by 
a similar hot spot [16]. By analogy with unwindings, the non-unwinding event is then taken 



to occur at the gridpoint around which the peaks are centred and at the time at which the 
cold to hot spot transition occurs. In every case the event occurs at the same gridpoint and 
on the same photon sheet. Since the sheets are separated by 1 grid spacing this specifies the 
event's time again to within 0.5 grid time units. 



4.3 CMB Anisotropies 

Finally we compare the CMB anisotropies induced by the texture field's evolution in each 
simulation. The details of the method are described elsewhere [fl6| 1; in essence we send a sheet 
of photons across a lattice over which the texture field is evolving and calculate the 'kick' 
given to each photon at each time step in the stiff source, small angle approximation [17]. 
The variation in the CMB temperature across the sheet is then recovered by Fourier analysis. 
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Here we are interested in the difference in the temperature anisotropics induced by a given 
initial field configuration evolved in each approach; we calculate 



1 N 

6T =nT, ((AT/T)f - (AT/T)? ) (23) 

i=l 

where the sum is over a given sheet of photons. For all configurations and for every sheet 
ST = 0. However the standard deviation of the temperature difference across the sheets 
clearly differs between the unwinding and non-unwinding configurations. Taking e = 8tt 2 G^, 
the magnitude of the anisotropy generated by the one known analytic solution to the flat space 
NLSM equations ||], we find 

f 0.05 e unwinding , . 

— 1 0.01 e non- unwinding 

Since the characteristic CMB anisotropy generated by an unwinding or non-unwinding 
event involves a sharply peaked cold spot being succeeded by a similar hot spot we can also 
consider differences in the maximum and minimum anisotropics between simulation pairs. 
We calculate 

1 N 

ST max = - £ ((AT/T)* ax - (AT/T)L X ) (25) 



1 N 



and 

5T miQ = -± £ ((AT/T) min - (AT/T) min ) (26) 
i=l 

where the sum is now taken over the set of all unwinding or all non-unwinding simulations, 
and now find 

g ( (-0.02 ± 0.04) e unwinding 

max | (—0.03 ± 0.03) e non-unwinding 

f (-0.04 ± 0.04) e unwinding 
min \ (-0.02 ± 0.02) e non-unwinding 1 ' 

where the error bars are la. Thus it appears that on average the NLSM gives marginally 
stronger peaks than the XCORE simulations, although in all cases the difference is within la. 
Indeed the only consistent discriminant between the anisotropies generated in each simulation 
is that the difference in individual peak heights (AT/T) max — (AT/T) m i n ) is in every case 
slightly greater in the NLSM simulation than its XCORE counterpart, although never by 
more than 0.1 e. 



5 Conclusions 

We have detailed the implementation of two approaches to circumventing the problem of 
the range of scales inherent in numerical simulations of global defect fields. It is to be 
expected that each approach will be constrained by a minimum lattice resolution, and that 



9 



such constraints will be strongest whenever significant events occur on microphysical length 
scales. In the case of global texture such events are the unwindings and non-unwindings 
associated with the local concentration of field gradient energy. We find that there is indeed 
a minimum resolution — of the order of 5 grid spacings per initial half-manifold radius - 
for such events to be modelled accurately for spherically symmetric simulations in flat space. 

For non-symmetric field configurations we have been unable to improve on the previ- 
ously established resolution limits [15] because of the difficulty of unambiguously identifying 
unwinding events in NLSM simulations of non-symmetric field configurations. The naive ap- 
proach to this problem is shown to be flawed in several ways, although the consequences for 
work in which it is used ju], [ll]] are unclear. It would certainly be possible to examine the 
effect of adopting inappropriate unwinding criteria on our single-texture simulations. How- 
ever because of the fundamantal differences — particularly in the grid resolution — between 
these and the many-texture simulations in which such unwinding criteria are used it would 
be difficult to meanigfully extrapolate any effect from one to the other. 

Working subject to these constraints we have been able directly to compare the evolution 
of a set of initial field configurations in each approach. In every case the fields are very strongly 
correlated, with what little discrepancy there is being almost entirely radial and simply 
reflecting the ability of the XCORE field to leave the vacuum manifold. Moreover, even 
where such discrepancies occur the subsequent return to near-perfect correlation indicates 
that the simulations are stable to them. Given this consistency it is not surprising that, 
insofar as we are able to identify them and to the accuracy of the identification schemes, all 
the unwinding and non-unwinding events occur at the same grid places and times in each 
approach. The fact that the fields remain so well correlated throughout the simulations also 
implies that the temporal derivatives of the field are equally consistent. Since the evolution 
of the velocity field that generates the CMB anisotopies their consistency — both absolutely, 
over individual simulations, and even more so statistically, over sets of simulations — is 
likewise a reflection of the correlation of the velocity fields. 

In conclusion we find that, subject to using a sufficiently fine lattice and a careful handling 
of the discontinuity inherent in the NLSM approach at unwinding, the XCORE and NLSM 
approaches to the simulation of the evolution of global texture are entirely consistent with 
one another, and hence may both be taken to give good representations of the true field 
evolution. 
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Figure Captions 



Figure 1 

a) Working on a 2-dimensional lattice consider each gridpoint i (empty dot) and its 4 nearest 
neighbours (filled dots) at time step n. 

b) Mapping the gridpoint and its neighbours onto the manifold the bisecting boundary is 
constructed by connecting the neighbours along geodesies (dotted line) . An unwinding should 
now be re-introduced whenever the field at the given gridpoint $j jn lies in the larger sub- 
section of the thus bisected manifold. The approximation considered here is now to replace 
the actual boundary with the 'average' boundary given by the intersection of the manifold 
with the plane perpendicular to the mean value of the neighbours <l>j jn (heavy line). Each 
boundary bisects the manifold into two distinct sub-sections. 



Figure 2 

For sufficiently (for example spherically) symmetric configurations $i >n is always parallel to 
<&i, n , so $i >n • $i jn < if and only if <&i n lies in the larger manifold sub-section — a) no 
unwinding, b) unwinding. In this case $j jn and $j jn being in different hemispheres is exactly 
equivalent to $j in lying in the larger sub-section. 

Figure 3 

In general <&j in is not parallel to 3>j in , and the equivalence noted in the symmetric case breaks 
down. Now can be in the same hemisphere as <&j in whilst remaining in the larger manifold 
sub-section. However replacing 3>j„ with $j >n — 3>j„ restores the equivalence, so the corrected 
unwinding criterion becomes ($i, n — ®i,n) ' &i,n < 0- 



Figure 4 

The approximation of replacing the true boundary with the 'average' boundary breaks down 
whenever $j n lies between the two. 

(a) if &i tn is in the smaller true sub-section but the larger 'average' sub-section then an 
unwinding will be mistakenly included. 

(b) if &i jn is in the larger true sub-section but the smaller 'average' sub-section then an 
unwinding will be missed. 



Figure 5 

The unwinding of a point other than the exact unwinding site in this prescription, where 
the field is unwound by (to first order) being inverted, 3>i )n +i ~ is clearly incorrect. 

For example there are configurations in which, if the neighbouring points do not evolve 
significantly, the unwound field remains in the larger manifold sub-section, and so re-unwinds. 
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Figure 6 

The variation in the topological charge Q causally associated with the unwinding of a spheri- 
cally symmetric texture field configuration with the grid resolution parameter TV; in both the 
XCORE (solid line) and NLSM (solid line) approximations. 



Figure 7 

The variation in the number of gridpoints less than 90% correlated between the XCORE and 
NLSM simulations of a typical random field configuration with run time. Both the full fields 
(solid line) and the normalised fields (dashed line) are compared. An unwinding event occurs 
in at grid time 13.9. 
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